Electronic structure of multiquantum giant vortex states in mesoscopic 

superconducting disks 
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We report self-consistent calculations of the microscopic electronic structure of the so-called giant 
vortex states. These novel multiquantum vortex states, detected by recent magnetization measure- 
ments on submicron disks, are qualitatively different from the Abrikosov vortices in the bulk. We find 
that, in addition to multiple branches of bound states in the core region, the local tunneling density 
of states exhibits Tomasch oscillations due to the single-particle interference arising from quantum 
confinement. These features should be directly observable by scanning tunneling spectroscopy. 
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I. INTRODUCTION 

Superconducting vortices are topological singularities 
in the order parameter |0| . In a bulk system each vortex 
carries a single flux quantum, while vortices with multiple 
flux quanta are not favorable energetically 1^. In small 
superconductors, however, the situation may be different. 
Today's nanotechnology can provide valuable insight into 
the nature of mesoscopic superconductors, whose linear 
dimensions can be comparable to the coherence length or 
the inter-vortex distance of the Abrikosov lattice. The 
following question then arises naturally: does single- 
quantum vortex matter survive the limit of decreasing 
sample size? More than thirty years ago, Fink and Pres- 
son gave the intriguing answer "not always" in their 
pioneering work Q on a related system: a thin cylin- 
der in parallel field. They have shown it theoretically, 
within the framework of the phenomenological Ginzburg- 
Landau (GL) theory, and also provided experimental evi- 
dence [|j for the existence of an enormous superfluid eddy 
current on the surface of a thin cylinder. They called 
this state a giant vortex state. While the work of Fink 
and Presson was largely forgotten for the next several 
decades, it nevertheless anticipated the present excite- 
ment in the field of nanoscale superconductivity. With 
recent advances in the controlled fabrication and study 
of nanometer-scale superconductors, the concept of giant 
vortex has been brought back to focus by Moshchalkov 
and coworkers a few years ago |^]. Their experiments 
on mesoscopic squares and square rings have indeed re- 
vealed that small superconductors do not always favor 
many-vortex states reminiscent of the Abrikosov vortex 
lattice. The measured H-T phase boundaries of these 
small structures were explained in terms of giant vortex 
states in the GL picture Subsequent experiments on 
submicron disks have further shown the existence of 



giant vortex states inside the phase boundaries. Within 
the GL framework |]6|-|l0| some of the abrupt changes in 
the magnetization observed have been attributed, e.g., to 
the collapse of a multi-vortex state into a giant vortex, 
or to transitions among different giant vortex states. 

This new phase of vortex matter has a single vortex 
occupying the sample, carrying multiple fluxoid |p|, pdHT2[ 
quanta. Such a state has no immediate analogue in 
bulk systems, and would only be similar to vortex states 
predicted for artificially patterned structures Jl3| , p^ . 
Moshchalkov et al. have also suggested that giant 
vortex states can cause the peculiar paramagnetic Meiss- 
ner effect seen in granular and mesoscopic superconduc- 
tors 16 ;^], through the compression of the flux trapped 
in the sample. This effect has been seen experimentally 
in mesoscopic systems by Geim and coworkers | p^ . As 
it is apparent from these recent experiments, mesoscopic 
superconductors exhibit novel quantum phenomena that 
are not observable in bulk systems. Studying their unique 
properties is crucial not only for potential applications 
but also for better understanding of nanoscale supercon- 
ductivity. 

Despite the fact that the existence of giant vortex 
states has been indicated more than three decades ago, 
and that their counterpart in nanoscale superconductors 
has been under intense scrutiny in recent years, there has 
been no microscopic, self-consistent theoretical study of 
its electronic structure. In this paper we report the re- 
sults of such a microscopic and self-consistent study of 
multiquantum giant vortex states in s-wave supercon- 
ducting disks of submicron size,^ using the Bogoliubov-de 



*After the completion of this work, we became aware of the 
recent results of A. S. Mel'nikov and V. M. Vinokur [^, in 
which the possible application of mesoscopic disks as quantum 
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Gennes (BdG) formalism . The spectroscopic proper- 
ties we have obtained for such vortex states can be probed 
directly by scanning tunneling microscopy (STM). Al- 
though GL studies give a good qualitative picture for a 
wide range of parameters, quantitatively reliable results 
are limited to the range relatively close to the critical 
temperature and magnetic field. Furthermore, analyzing 
the system from a microscopic point of view is essen- 
tial to understanding superconductivity on such a small 
scale. Latest experimental efforts aimed at STM imag- 
ing of mesoscopic vortex matter have focused on NbSe2 
samples since direct STM images of vortex states 
have been obtained only on high-quality single crystals 
of NbSea (H] and BiaSrsCaCuOg+a Hll. These com- 
pounds are highly two-dimensional and easy to cleave 
in situ, providing very clean surfaces - a key ingredient 
for successful STM imaging. Thus, anticipating STM 
measurements, we present in this paper self-consistent 
BdG results with parameters corresponding to submicron 
NbSe2 disks. 



II. FORMULATION 

We consider an s-wave superconducting disk of radius 
R under a magnetic field perpendicular to the disk area, 
with a vortex carrying m fluxoid quanta formed in the 
center. The system has cylindrical symmetry and it is 
described using cylindrical coordinates [r,9,z). In ac- 
cordance with the experiments |^,^ we assume the disk 
thickness to be much smaller than the penetration depth. 
Consequently, the order parameter is assumed to be uni- 
form in the field direction (z), and the current density 
j as well as the vector potential A has no z component 
0. In the gauge which removes the phase of the order 
parameter i.e., A(r) = |A(r)|, we can write down 
the radial part of the BdG equations pMlSpq] as 
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+ (7,{U{r) - n)iJn{r) + a^\A{r)\iPnir) = e„V'«W, (1) 



where 



u„{r) 
v„{r) 



(2) 



is the radial quasiparticle amplitude. Here crx and are 
the Pauli matrices, rUe is the electron mass, and /i the 



switches is studied. This work discusses the tunneling density 
of states in the core of giant vortices in terms of quasiclassical 
calculations. 



chemical potential. The angular momentum quantum 
number I is an integer when m is even, and a half odd 
integer when m is odd . The single-particle potential 
U{r) can incorporate the lattice potential and inhomo- 
geneity effects due to impurities and sample boundaries. 
To study quantum size and interference effects, we con- 
sider a clean sample and take U{r) — inside the disk, 
while including the periodic lattice potential in terms of 
the effective masses, and m^. Furthermore, we take 
ruz 3> rur, as justified for highly anisotropic materials 
such as NbSe2 , and neglect the dependence of Eq. on 
the motion along the z direction. 

In principle, due to finite thickness of the sample, the 
vector potential must depend on not only r but also z, 
so that A = Ag{r, z)6, and Ag{r) in (Q) is an average of 
Ag{r, z) over the disk height L 



1 /-^/^ 
Mr) = T / dzAe{r,z). 



L 



(3) 



-L/2 



However, for typical experimental parameters the 
lateral sample size is much larger than the thickness, and 
we therefore consider the case where the vector potential 
is independent of the z coordinate; Ag(r, z) ~ Ao{r). The 
order parameter and the current density j = je{r)9 are 
given in terms of the eigenvalues and eigenfunctions of 
Eq. (0) as 



|A(r)|=g Mr)v:ir){l-2U) (4) 
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(5) 



where g is the coupling strength for the electron-electron 
attraction, ujo is the Debye frequency, and /„ = f{en) is 
the Fermi distribution function. The vector potential is 
given in turn by the current density through the Maxwell 
equation V x V x A = (47r/c)j. We first solve Eq. (|l|) 
with initial guesses for |A(r)| and Ag{r) and recalculate 
them from Eqs. (^) and (H), and repeat the process until 
self-consistency is acquired. The local tunneling density 
of states and the differential conductance [ p5|JT^ can then 
be calculated by 

7V(r,i?)-^[|zi„(r)|25(i?-e„) 

n 

+ \v^{r)\H{E + e^)] (6) 

n 

+ l^'n(r)P/'(e„ + ey)]. (7) 

Clearly, the differential tunneling conductance (|^) is a 
direct probe of the local density of states, N{r,E), pro- 
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vided that the temperature is low enough. In the fol- 
lowing we will present some results for the experimen- 
tally observable differential conductance at low temper- 
atures, but may refer to it as the local tunneling den- 
sity of states (LDOS). We choose the parameters corre- 
sponding to NbSe2: we take — 2me, Ep = 37.3 meV, 
hwD = 3.0 meV, and set the coupling strength so that 
the bulk gap Aq = 1.12 meV. 

The results are shown for disk radius R = 500 nm. We 
have investigated the size range of i? = 200 — 600 nm and 
have found qualitatively same features in the LDOS. 



III. RESULTS 
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FIG. 1. Local tunneling conductance as a function of coor- 
dinate r and voltage V, for a giant vortex state with (a) m = 4 
and (b) m = 5 flux quanta, sustained in a superconducting 
disk with radius R = 500 nm at temperature T = 1 K. 

In Fig. 1 we show the differential conductance for a 
vortex with (a) m = 4 and (b) to = 5 flux quanta, as 
a function of voltage V and radial distance r from the 



disk center, for temperature T = IK []. In both cases, 
prominent sharp peaks can be seen near the vortex core 
and for low voltages - four peaks in the former and five 
in the latter. Generally speaking, the number of low-bias 
conductance peaks corresponds to the winding number 
TO of the order parameter, which gives rise to to peaks 
near the center This feature is in accordance 

with the index theorem established by Volovik for the 
Caroli-de Gennes-Matricon bound states in a vortex core. 
According to this theorem, the quasiparticle spectrum 
of a vortex with winding number m has m branches of 
bound states, which cross zero energy as a function of 
angular momentum. These quasiparticle branches also 
explain the evolution of the m rows of conductance peaks 
as one moves away from the core, as seen in Fig. 1, with 
decreasing number of peaks one by one fl^ . In contrast 
to the singly quantized case |28|, one can see directly 
that as energy increases, more states with higher angular 
momenta contribute to the density of states. 

We illustrate this in Fig. 2, where a spatial map of the 
LDOS is taken for various fixed values of V. Clearly, 
with increasing bias voltage the density of states is redis- 
tributed from the core towards the sample boundaries. 
Fig. 1 also reveals the presence of the so-called zero mode, 
i.e., a peak around zero energy at the vortex core for 
odd TO, and its absence for even to. The existence of 
a zero mode is, quite generally, linked to a sign change 
in the order parameter as a function of some general- 
ized coordinate | p9| . The zero-bias peak is a signature 
of bound states for quasiparticles trapped by the sign 
change. Here, in the given gauge, the order parameter 
changes sign at the vortex core when m is odd, while it 
does not when m is even. 

The low-bias peaks and zero modes discussed above 
are general characteristics of the LDOS associated with 
the winding number of the order parameter. In addi- 
tion to these, however, we have found novel features in 
the LDOS that are unique to giant vortex states in sub- 
micron disks. They are the oscillations seen above the 
gap energy in Figs. 1 and 2(c) and more clearly in the 
contour plot of Fig. 3. These oscillations are similar in 
origin to the so-called Tomasch oscillations discovered in 
a superconductor-normal metal junction [^0|-^ and re- 
flect "standing waves" arising from the interference of 
quasiparticle states. This interference effect is a direct 
consequence of strong confinement experienced by the 
superconducting quasiparticles due to the small system 
size. We dedicate the remainder of this paper to detailed 
discussions of this effect. 



^For data to be shown for energies well above the gap energy, 
fast l/kp oscillations have been removed by means of Fourier 
transform. In actual observations these fast oscillations will 
not be resolved. 
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FIG. 2. Spatial map of the local tunneling conductance for 
the entire disk in the giant vortex state of Fig. 1(b) , for various 
fixed voltages. It can be seen that the maxima in the local 
density of states gradually shift towards the perimeter of the 
disk as the voltage is increased. 

When a vortex holds multiple flux quanta m, the or- 
der parameter vanishes around the center over a cer- 
tain area - the larger the m is, the larger the area [ p6[ . 
As the distance from the center increases, the order pa- 
rameter increases and recovers to its bulk value eventu- 
ally. In the case of NbSe2, due to the short coherence 
length, the recovery happens relatively quickly. This re- 
sults in a well-defined superconducting region with the 
constant order parameter Aq within the disk. At the 
disk boundaries, however, the order parameter is forced 
to vanish and as a result, exhibits Friedel-like oscilla- 
tions around its bulk value near the surfaces. These 
oscillations have the largest amplitudes at the bound- 
aries and decay roughly over the coherence length scale. 
Moreover, the smaller the system size, the larger these 
amplitudes relative to the bulk value. The quasiparti- 
cles confined in the disk experience scattering by this 
large change in the order parameter at the surfaces. An 
electron-like quasiparticle is reflected back as a hole-like 
one and vice versa, and the Tomasch effect results from 
the interference between the electron-like and hole-like 
states in the superconducting region js^]. The momenta 
of electron-like and hole-like quasiparticles for energy E 
are = ^^^^^/Ep ± ^ kp ± f§^, respectively, 
where fl = E"^ - |A(r)|2 with |A(r)| ~ Aq and kp is 
the Fermi momentum. At a given distance d from the 
surface, the LDOS oscillations in energy are determined 

by H 1^ ^ + , where n is an in- 

teger, and ~ 151.15 nm for NbSe2. Furthermore, 

the interference can be seen in the LDOS also as a func- 
tion of coordinate (distance from the surface) for a given 
energy E. The period of the oscillations in this case is 
given by 8d ~ i^fivF i ^j^g oscillation peri- 

ods obtained in our numerical results, as seen in Fig. 3, 



are in quantitative agreement with these analytical ex- 
pectations. 




FIG. 4. Same as Fig. y, but for the corresponding model 
calculation. 

In a superconductor with short coherence length, if 
the winding number m and, consequently, the "normal 
core" area is very large, the Tomasch effect may arise 
also from the vortex core, as in a normal-superconductor 
junction. For submicron NbSe2 disks with m up to five, 
however, we have found that the LDOS is dominated by 
the Tomasch oscillations coming from the surfaces. In- 
deed, we have confirmed the LDOS oscillations charac- 
teristic of the Tomasch effect in terms of model calcula- 
tions in one and two dimensions, where the BdG equa- 
tions are solved without iteration with a step-function 
order parameter: A(r) = 0(r < i?/2);Ao(r > i?/2). In 
this case the Tomasch effect coming from the normal- 
superconductor interface governs the LDOS structure, so 
that the oscillation period in energy (see En above) be- 
comes larger as one approaches the interface (i.e., here 
d is the distance from the interface). Apart from this 
and enhanced amplitudes due to a larger change in the 



4 



order parameter, the LDOS shows the same quahtative 
features as seen above (compare Figs. 3 and 4). 

IV. CONCLUSIONS 

We have presented detailed, self-consistent calculations 
of the microscopic electronic structure of giant vortex 
states. We believe that the most direct experimental ev- 
idence for the existence of giant vortices can be provided 
by STM measurements of the local density of states in 
sub-micron superconductors capable of sustaining such 
vortex configurations. We have provided a spatial map 
of the tunneling density of states for multi-quantum gi- 
ant vortex states, and have identified several signatures 
that can be used to identify them with STM. We have 
found that under extreme confinement the quantum in- 
terference arises among quasiparticle states and leads to 
experimentally observable Tomasch oscillations in the lo- 
cal density of states. 
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